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Abstract 

The motion of massive test particles in dark matter is studied. It is shown 
that if the energy density of the dark matter making up a galactic halo has a 
large r behavior of then contrary to intuition the motion of these test 
particles are not govern by Newtonian gravity, but rather by the equations 
of geodesic motion from Einstein's theory of general relativity. Moreover, 
the rotational velocity curves of orbiting massive test particles in this energy 
density do not approach a constant value at large r but will instead always 
increase with the radius of the orbit Tc- 
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§1. Introduction 

One of outstanding problems in astrophysics today is to explain the 
motion of bodies with orbits in the galactic halo. Prom Newtonian dynamics 
one would expect that for a body in a circular orbit 

^ GMg 

where G is the gravitational constant, Mg is the luminous mass of the galaxy, 
Tc is the radius of the circular orbit and is its rotational velocity. If one 
were then to plot verses r for various bodies orbiting in a galactic halo 
one would expect to obtain a rotational velocity curve which decreases as 
l/V^c- ^^^^ does not, in fact, happen. Experimentally, it is instead found 
that the rotational velocity curve approaches a constant value at large r. To 
explain this result, the presence of "dark matter", matter which has yet to 
be detected, has traditionally been proposed. (See [1] — [3] and the references 
contained therein.) Namely, it has been postulated that the galactic halo is 
filled with a gas of weekly interacting particles with a mass density p. With 
their presence the total mass contained within the radius of the orbiting 
body changes and the Newtonian equation of motion now becomes 
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If we now take p ^ aoc^/{Gr'^) for large r where ao is a dimensionless 
constant, then as long as 

MgG 

the Mg term may be neglected and a constant f <^ with value 

vl 

-| 47rao , 



can be obtained. (Rg is the point at which the velocity curves become a 
constant and is identified as the radius of the luminous galaxy.) Explana- 
tion of the rotational velocity curves then reduces to finding and detecting 
candidates for this dark matter which will have the correct large r behavior. 

The basic premise of this argument is that Newtonian dynamics and 
gravity will still be valid even after the introduction of a energy den- 
sity into the system. This need not be true. One should remember that 
Newtonian gravity is only an approximation of Einstein's theory of general 
relativity and that the very act of introducing a mass density p ~ for 
dark matter introduces an unconfined energy density into the system. Its 
presence cannot help but have an affect on the geometry of the spacetime 
in the halo, and, consequently, on the motion of bodies orbiting there. 

In this paper we shall show that the standard argument using New- 
tonian dynamics and gravity for the existence of dark matter with a 
energy density is inconsistent. The introduction of an energy density which 
behaves as 1 /r^ for large r changes the geometry of the spacetime so dras- 
tically that the motion of bodies in the galactic halo is necessarily non- 
Newtonian. The basic premise that Newtonian gravity is still valid even 
after the introduction of this dark matter is incorrect. Furthermore, after 
using the full geodesic equation from Einstein's theory of general relativity 
to analyze circular orbits in the energy density, we find that v'^ ~ qr^ 
for < q < 1. The rotational velocity always increases with Tc and will not 
approach a constant value. In fact, for a static, spherical geometry there 
are only a very narrow range of essentially unphysical choices for the energy 
density which will give a constant for r > Rg. 

§2. Geometry of the ao/r^ energy density 

We begin by modeling the galaxy as a sphere of mass Mg and radius 
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Rg which is surrounded by a galactic halo made up of dark matter. The 
most general static, spherically symmetric metric is known to be [4] 

ds^ = -f{r)dt^ + h{r)dr^ + r^{dd^ + sin^ Odcj)'^) , 

where / and h are unknown functions of r only which need to be determined. 
As usual, we write the energy momentum tensor for the system as 

where p is the energy density of the particles in the halo, p is their pressure, 
and Ufj^ is an unit velocity vector in the direction of the timelike Killing 
vector for the system. From Einstein's equations [4] we then obtain 

87rp- — ^ r A _ 1 
^ rhf ^ h 

« ir 1/' ff\h'\ f h' 



2fh 4fh\f hJ 2rfh 2rh? ' 
where we are using units in which G = c = 1 and the primes denote deriva- 
tives with respect to r. 

Suppose now that the energy density of dark matter is p ^ ao/r'^ for 
r > Rg where ao > is a (dimensionless) constant. Then the first equation 
in (2) is trivial to integrate giving 

, _i K 

h — 1 — 8nao , 

r 

where K is an integration constant. Since in the absence of dark matter 
we would expect to have obtained the Schwarzchild solution, we identify K 
with 2Mg. Next, notice that as long as the Newtonian bound (1) on ao 
holds 

_^ Ma 



and we can neglect this term in h and can approximate ~ 1 — 87rao as 
a constant for r > Rg. We shall justify this approximation later. 
The difference of the second two equations gives 



2/ 4f\f^hJ 2r\f^hJ^ r2 



Since h is a constant for r > it is straightforward to solve, 

fir) ^ (k+r'i+/^ + k_r'i-/^y , 

where 



g± = 2 ± 2 



V 1 - 87rao ) 



and k± are integration constants. Because / > 0, ao < l/167r otherwise 
/ will contain oscillatory solutions. Since we are working in the large r 
limit, and since q+ > q^, only one of the two solutions to (2) will survive. 
Linear combinations of the two will not. We can thus consider each solution 
independently and for convenience we shall write f± = k±r'^^. Then 

q± + l-h 

and we see that the pressure also varies as Moreover, 



p 



4-g± 



and since > 2, p+ > p. Because p < p/3, we must therefore exclude 
the (+) solutions as being unphysical. The only physical solutions are the 
(— ) solutions, and we thus set q = q^, p = and / = kr'^. Furthermore, 
< g < 1 while Mg/{A7rRg) < ao < 3/(567r). Consequently, 1 < h < 7/4, 
and, since / ~ r^, we can see explicitly that spacetime in the presence of 
the energy density is not flat, but is instead quite curved. 
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We next consider the motion of a massive test particle in a static, 
spherically symmetric geometry with a velocity such that — 1 = Vf^v^. 
This constraint gives 

where r is the proper time of the particle. Working in the equatorial ^ = 7r/2 
plane, and using energy and angular momentum conservation, 

where E and L are the energy and orbital angular momentum per unit mass, 
respectively, of the particle, we obtain 

2 

dt 



0=(?) +V{r), 



where 



is an effective potential energy. For circular motion, r = Tc, the radius of 
the circular orbit which is a constant in time. Consequently, V{rc) = and 
V'{rc) = 0, giving 

2/(re) - re/' (re) ' r2 2f{r,) - rj'{r,) ' ^ ^ 
Defining the rotational velocity as 

d(f) 

then from (4) and (5) we find that for circular motion. 



vl = ^re/'(re) . (6) 
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This equation holds for any /. In particular, if we are dealing with a spher- 
ical mass M in free space, then / = (1 — 2M/r) and (6) reduces to what 
one obtains from Newtonian gravity. If, on the other hand, p ~ 1/^^, then 
from the above / = /cr^, so that 



For q ^ 0, Vff) always increases with the radius of the orbit and never ap- 
proaches a constant value as one would naively expect from Newtonian 
gravity. 

In the above solution of (3) for / we have neglected the contribution 
of h'/h in comparison to /'//• Wc shall now justify this approximation. 
First, we note that for the galaxy, 2M/r <^ 1 and, due to the bound (1) 
on ao, this term is very small in comparison to 1 — Sttoo. Next, note that 
while h'/h ~ 2M/r^, the solution we obtained by taking /i as a constant 
gives /'// ~ I/'''- Since r is large, we would expect h'/h to have a very 
small affect on the kr"^ solution of (3). Consequently, we can take h'/h as 
a small perturbation and solve (3) perturbatively about the kr'^ solution. 
After doing so, we find that to first order in 2M/r, 



where hg^ = 1 — Sttoo. The inclusion of the 2M/r term in h modifies the kr'^ 
solution very slightly since 2M/r <^ 1, and 1 < /lo < 7/4. Consequently, we 
were justified in neglecting this contribution to /. 
§3. Constant v^f, energy density 

We now ask whether or not it is possible for any physically reasonable 
p to result in a which will be constant for r > Rg. Using (6), we find 
that for V(j, is to be a constant outside of the galaxy, / must then have the 



v. 
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approximate form of fy 2v|log(r/ro) for r greater than some vq. (The 
subscript v will denote the fact that we are looking for solutions of Einstein's 
equations which will result in a constant f^.) Because > 0, r > tq, and 
since is a constant only outside of the galaxy, we shall identify tq with Rg. 
Since /„ is now given and unknown, (3) becomes a differential equation 
for hy which may be written as 



dy \hyj y \hy J 1 + 2y ' 
where y = log(r/i?g) > 0. Its solution is straightforward, 

where is an integration constant. Then from (2) we find that: 

=1 + ^ - (1 + ,y)e^y (c. + 2 £ (l-^dt) , 
8nr^py^{l + y)e^y (^Ch + 2j^ ^^d?j-l. (7) 

Physically, p„ > 3p„. This gives an upper bound of 0.0013 for c/j. As /t > 
for all J/ > 0, c/i > 0. Consequently, < Ch < 0.0013 and there are only a 
very narrow range of values for Ch which will result in a physically reasonable 
Py and py. Moreover, if > 0, then it is only when y is between some ymin 
and ymax that py > 3py. Outside of these two bounds the energy density 
must have a different form and cannot be a constant. If, on the other 
hand, c^, = then py > 3py for all y > 0.627. As long as r > e^-^^^i?^ the 
energy density given in (7) for Ch = will result in a constant rotational 
velocity curve outside of the galaxy. For r < e^-^'^'^Rg the energy density 
will have a different form and V(j, will not be a constant, as expected. 
§4. Concluding Remarks 



We have thus shown that if the energy density p of dark matter behaves 
as for large r, then / fa fcr^ for some constant k while h 1 — Sttoq. 
Consequently, the premise that p ~ for large r contradicts the premise 
that Newtonian gravity is valid in the halo. The introduction of this energy 
density, which is not confined but spread over a large area, alters the geom- 
etry of spacetime so drastically that Newtonian dynamics is no longer valid 
in the halo. The assumption that Newtonian gravity is valid even after the 
introduction of dark matter with a energy density is incorrect. In fact, 
contrary to what is expected from Newtonian gravity, circular orbits for this 
energy density have a v| = qkr1/2 which always increases with rc- It never 
approaches a constant, although for very small q it increases very slowly. 

From a physicist point of view, however, the major problem with using a 
energy density for dark matter is not that there was an inconsistency in 
the argument for its introduction since this is straightforwardly resolved by 
using general relativity instead of Newtonian gravity to analyze the system; 
nor is it that f ^ always increases with Tc since q may be taken to be quite 
small so that any increase in Vcf, occurs very gradually (with one cavet; 
see [5]). It is rather that the spacetime in a energy density is so 

curved that Newtonian dynamics is no longer valid. If the energy density of 
the dark matter in the galactic halo truly does behave as for large r, 
this will present great difficulties in interpreting extragalactic astronomical 
observations. As most galaxies have a halo, including presumeably our 
own, light from a distant galaxy would then first have to pass through its 
own halo, a region of curved spacetime, and then our through own galactic 
halo, another region of curved spacetime, before we can observe it. Since 
the usual assumption is that light from other galaxies passes through a 
spacetime which is essentially flat before it reaches us, if the energy density 
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of dark matter has a large r behavior, then all of the extragalactic 
observational data would have to be re-evaluated and interpreted. This is 
just one of the many problems that would arise from using a energy 
density for dark matter which have yet to be addressed. 

The energy density is an unconfined energy density which, pre- 

sumeably, extends for large distances into the spacetime. As the universe 
is known to be expanding and thus changing with time, one may question 
the validity of using a static solution of Einstein's equations to analyze the 
motion of test particles in the dark matter as we have done. As, however, 
the galactic rotation curves only extend out to a few galactic radii, we are 
only interested in the behavior of the motion of bodies relatively close to 
the galaxy and in this region a static approximation is certainly valid. Of 
course, at some r ^ Rg the 1 /r^ energy density must be cutoff-ed and our 
analysis will no longer be valid, but this will only happen outside the region 
we are interested in. We should also note that precisely the same static ap- 
proximation is made in the standard analysis of the motion test particles in 
the galactic halo using Newtonian gravity. The important point here is not 
whether the use of a static, energy density to model the galactic halo 
is a valid approximation or not, but rather that our analysis of the motion 
of test particles in a energy density using general relativity holds in 
precisely the same regime in which the standard analysis using Newtonian 
gravity was presumed to be valid. 

As (3) is a non-linear second order differential equation in /, it may 
be that we truly cannot neglect the h' /h term in (3) no matter how small 
2M/r is. In §3, however, we have shown that in a static, spherical geometry 
a rotational velocity curve which is truly a constant for r > Rg is obtainable 
only for the very special choices of the py given in (7). The only approx- 

10 



imation made in obtaining (7) was once again that the system is static, 
and spherically symmetric. Notice, however, that if (7) is truly the energy 
density of dark matter, this would mean that no matter what the internal 
properties of the galaxy are, once one leaves it the energy density of the 
particles making up its' halo has to be determined within one part in 1000. 
This is extraordinarily and prohibitively restrictive. The form that py takes 
is very particular and it is difficult to imagine a physical process which will 
not only reproduce it, but also determine py to such a high degree of ac- 
curacy. On this basis alone we would tend to rule out py as a physically 
viable energy density for dark matter. We also note, however, that because 
/ ~ log (r/i?g), once again the spacetime with this energy density is curved 
and in using py as the energy density of dark matter we would once again 
be faced with the problem of interpreting the observational data. 

The analysis done in this paper was done for a very special system un- 
der some very restrictive conditions. For example, although we have used a 
spherically symmetric geometry to model the galaxy, most observed galaxies 
are axisymmetric and have a definite angular velocity. It is, moreover, not 
even clear whether the experimental data gives rotational velocity curves 
which are truly a constant or whether they are instead slightly increasing or 
decreasing with rg. All that we are comfortable concluding from this analy- 
sis, therefore, is that one must be much more careful about introducing any 
unconfined energy density for dark matter into the system. Not only must 
it be able to explain the experimental rotational velocity curves within the 
framework of general relativity, but one must also consider the subsequent 
affects of this energy density on the geometry of the spacetime. As we have 
seen, for the 1/r^ energy density these affects are considerable. 
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